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High precision measurements of the renormalized zero-momentum 4-point coupling ga and of the Liischer-Weisz- 
Wolff running coupling g(L) — Lm(L) performed with two different lattice actions in the non-perturbative region 
confirm the earlier observations, that the cutoff effects look linear, in contrast to perturbative considerations. The 
use of different actions allows one to make a more reliable estimate on the continuum limit. The measurements 
were done for infinite volume correlation length up to 350. 



In numerical simulations with lattice regular- 
ization the physical results are obtained after ex- 
trapolating to zero lattice spacing. Usually (and 
especially in 4d) the range of a values accessi- 
ble in the simulations is quite restricted, there- 
fore the extrapolation a — > is highly non-trivial 
and the knowledge of the functional form of the 
artifacts is essential. As shown by Symanzik |Q, 
in bosonic theories in every order of perturbation 
theory (PT) the leading lattice artifacts decrease 
like 0(a 2 ) (apart from logo factors). It is gen- 
erally assumed that this behaviour holds beyond 
PT, and the extrapolations in numerical simula- 
tions are done according to this assumption. 

In Q the renormalized zero-momentum 4-point 
coupling was measured in the 0(3) non- linear 
sigma model to high precision. Unexpectedly, the 
cutoff effects were nearly linear in a and were de- 
scribed well by a form cq + c\a + c 2 a log a, in con- 
trast to the standard belief. (Note that already 
in Q there was some numerical evidence that a 
linear Ansatz describes the data better, but the 
errors were significantly larger.) 

Here we extend the results of [|| to another 
quantity, the Liischer-Weisz- Wolff running cou- 
pling g(L) We also measure these quantities 
using an alternative (ad hoc) lattice action con- 
taining a diagonal interaction term. For this ac- 
tion the coefficient of the nearest neighbour (on- 
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axis) term is (3\ while the coefficient of the diag- 
onal term is fti- They were chosen to be equal, 
&x = 02 = (In this notation for the standard 
action fi\ = f3, fii = 0.) 

The simulations were done with both actions at 
the same values of the correlation length and lat- 
tice size (i.e. same a/L and £(£))■ The assump- 
tion of universality implies that both actions lead 
to the same continuum limit for physical quanti- 
ties. This strongly restricts the fits to the data. 

The quantities measured are: 

1) The renormalized 4-point coupling at finite 
volume 

/ L \ (5 ((M 2 ) 2 ) \ 

defined on an L x L periodic lattice. Here M = 
J2 X S is the total magnetization while ^(L) is the 
second- moment correlation length (cf. e.g. ||). 
The physical size, z = L/£ 2 (L) is kept fixed. 

2) The running of the LWW coupling (§. The 
coupling g(L) = m{L)L is defined on an L x oo 
strip, where m(L) is the inverse correlation length 
in the zero-momentum spin-spin correlation func- 
tion. For a given number L/a of lattice sites in the 
spatial direction one tunes the value of (3 to get a 
fixed value g(L) — uo. In the next step one mea- 
sures g(2L) at the same value of (3, giving g(2L) = 
£(2, it , a/L) The physical (regularization in- 
dependent) "step scaling function" is given by the 
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gn(zo) for the two actions at zo 
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3.05 
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(mod) 


5 
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1.33591(12) 
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1.27668(9) 


1.29060(13) 
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1.27228(12) 


1.28242(21) 
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1.26817(9) 


1.27424(10) 
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1.26306(16) 


1.26456(20) 


Table 2 






£(2,u 


a/L) at i*o = 


1.0595. 



limit lim a ^o S(2, uq, a/L) = a(2,uo). Due to an 
improved estimator by Hasenbusch || (suited for 
the strip geometry) we could significantly improve 
the error. 

The values of 3_r(zo) a t z q — 2-32 for different 
lattices sizes L/a and the two actions are given in 
Table [l]. These data are shown in Figure |l|. Ta- 
bic H and Figure refer to the LWW step scaling 
function Y<(2,uo,a/ L) at uo = 1.0595. 

The data with L/a > 10 were fitted simultane- 
ously for the two actions by several forms: 

a a L 
A: c + ci— + c 2 — log— , 
L L a 

B: C o +Cl (f) 2 + c 2 (f) 2 log^ (2) 

The coefficients Co, c±, ci correspond to the 
standard action while those for the modified ac- 
tion we denote by c , c' 1i c' 2 . Universality implies 
c' = Co, and for the case C we also took the same 
power, c' 2 = C2- Table || gives the coefficients and 
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Figure 1. gji(z) at z = 2.32. The data for the 
larger artifacts correspond to the modified action. 
The fit shown has a and a log a terms. 



the values of \ 2 /dof for the fits. (The pure a 2 fits 
have x 2 /dof = 11800/7 and 1270/7 hence are not 
listed.) 

Interestingly, one does not necessarily need to 
assume a specific analytic form of the cutoff de- 
pendence. Since we performed the measurements 
with the two actions at the same values of a/L 
and it is enough to assume that the form 

of the leading cut-off effect is the same, only the 
prcfactors are different. A quantity q measured 
with the two actions is given then by 

q(a) = q + f(a) + ... (3) 
q'{a) =q + cf(a) + ... 

where /(a) is an unspecified function. We assume 
that at the values of a considered the non-leading 
artifacts are negligible. Measuring the quantity 
q at different values ai,et2,... , a^r one has 2N 
data q(a.i), q'(di) for the two actions and only 
A^+2 unknown variables in the fit, qo, c and fifli). 
Both gn(zo) and cr(2,uo) are described well by 
the assumptions made in (|^) . In Table |] we give 
the continuum value, c and x 2 /dof for the two 
quantities measured.^] 

For small g(L) the cutoff effects can be stud- 
ied in PT. The sign of artifacts at g(L) = 0(1) 

2 Obviously, by taking three different actions one can re- 
solve the next-to-leading cutoff term as well. 
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Figure 2. The step scaling function S(2,Mo,a/L) 
at Mo = 1.0595. The data for the larger artifacts 
correspond to the modified action. The fit con- 
tains a and a log a terms. 



(like in Fig. g) is opposite to that obtained in 
PT. This has been noticed already in where 
it has also been pointed out that the large N 
limit shows the same effect. Caracciolo et al. 
(g| studied the large N integrals analytically 
and found (a/L) 2 (log L/a) q artifacts with q = 
— 1, 0, 1, 2, . . . in the step scaling function. They 
pointed out that the negative powers of log L/a 
are not reproduced by PT. 

Based on our numerical study in 0(3) it can not 
be ruled out that for very small a the perturbative 
prediction (a 2 times powers of log a) takes over. 
In this case an explanation would be needed why 
this happens at £ > 350 only. 
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Table 3 
Fits of the form 



for g R (z ) and cr(2, u ) 





9o 


c 


X 2 /D 
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Table 4 

Results of the form-independent fits. 
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